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ABSTRACT 

The ejection of gas out of the disk in late-type galaxies is related to star formation 
and is mainly due to the explosion of Type II supernovae. In a previous paper, we 
considered the evolution of a single Galactic fountain, that is, a fountain powered by 
a single supernova cluster. Using 3D hydrodynamical simulations, we studied in detail 
the fountain flow and its dependence with several factors, such as the Galactic rota- 
tion, the distance to the Galactic center, and the presence of a hot gaseous halo. As a 
natural follow-up, the present paper investigates the dynamical evolution of multiple 
generations of fountains generated by ~ 100 OB associations. We have considered the 
observed size-frequency distribution of young stellar clusters within the Galaxy in or- 
der to appropriately fuel the multiple fountains in our simulations. Most of the results 
of the previous paper have been confirmed, like for example the formation of interme- 
diate velocity clouds above the disk by the multiple fountains. Also, the present work 
confirms the localized nature of the fountain flows: the freshly ejected metals tend to 
fall back close to the same Galactocentric region where they are delivered. Therefore, 
the fountains do not change significantly the radial profile of the disk chemical abun- 
dance. The multiple fountains simulations also allowed to consistently calculate the 
feedback of the star formation on the halo gas. We found that the hot gas gains about 
10 % of all the SNII energy produced in the disk. Thus, the SN feedback more than 
compensate for the halo radiative losses and allow for a quasi steady-state disk-halo 
circulation to exist. Finally, we have also considered the possibility of mass infall from 
the intergalactic medium and its interaction with the clouds that are formed by the 
fountains. Though our simulations are not suitable to reproduce the slow rotational 
pattern that is typically observed in the halos around the disk galaxies, they indicate 
that the presence of an external gas infall may help to slow down the rotation of the 
gas in the clouds and thus the amount of angular momentum that they transfer to the 
coronal gas, as previously suggested in the literature. 
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1 INTRODUCTION 

Galactic fountains (GFs) are believed to occur in the Milky 
Way as well as in the other disk galaxies. A GF takes 
place when the Type II supernovae (SNe II) belonging to 
an OB association are sufficiently numerous to create a su- 
perbubble and to dri ve its expansion above the scale height 



of the gase ous disk (|Shapiro fc Fieldl 1 1976t Bregman 
lKahnlll98ll ). As the break out occurs (Ma c Low et al 
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iKoo fc McKeell 19921 ). the gas heated by the SNe II falls back 
to the Galactic plane as cold gas, mostly concentrated in 
clouds formed through thermal instabilities. 

The GF mechanism is thought to be linked to a num- 
ber of issues which have been discussed in some detail in 
iMelioli et~ al, (2008, thereafter Paper I). Here we briefly re- 
call them: i) the formation of the observed thick H I layer 
having a mass of 1/10 of the total H I, and rotating with 
velocity about 20-50 km s _1 smaller than that of the gas on 
the plane; ii) the presence of high velocity clouds (HVCs) 
and/or intermediate velocity clouds (IVCs). It is important 
to understand whether these clouds are forme d by GFs, or 
represent external gas accreting on the disk (jSancisi et al.1 
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2008); Hi) the presence of giant holes in the H I disk of 
many galaxies, which can be produced by associations of 
O and B stars exploding in a correlated fashion and giv- 
ing rise to superbubbles; iv) the chemical evolution of the 
Galaxy, as in principle GFs might be able to move SN II 
ejecta relatively far from the place where they have been 
produced, affecting the gradients of a-elements; v) the dif- 
fuse, soft (T < 1 keV) , extraplanar X-ray emission detected 
in some spiral galaxies which is due to hot gas extending up 
to 10 kpc above the galactic disk. This gas is most likely 
blown out by star forming regions in the disk, as the X-ray 
luminosity correlates with the star formation rate and the 
hot gas is oxygen rich, a signature of SN II enrichment. 

Given the richness of the phenomena involved in the 
GFs, a number of papers has been devoted to this argu- 
ment. A detailed analytical description of the motion of a 
fluid element in a GF from the rotating disk to the halo and 
back, predicting the locati on where it w ould return on the 
disk was first presented by iKahn lll99ll) and Kahn fc Bretd 
(|l993h . lde Avilled (|2000h and lde Avillez fc Berrvl (|200ll ) per- 
formed several 3D hydrodynamical simulations of the gas in 
the Milky Way in order to account for the collective effects 
of supernovae on the structure o f the interstellar medium 
(ISM). In particular, Ide Avilled (2000) discussed in detail 
for the first time the fraction of material that results in the 
formation of IVCs. Mor e recently, further 3D MHD simula - 
tions were carried out bv lde Avillez fc BreitschwerdtJ (|2005h . 
They found that gas transport into the halo as well as the 
mean volume filling factor of the hot phase in the disk is 
not significantly affected by t he presence of a initial disk 
parallel magnetic field (see also iTom isaka 1998I). In order to 
reach a high spatial resolution, all these simulations consider 
only small disk regions (typically 1 kpc 2 x 10 kpc) an d as- 
sume plane parallel stratified ISM. iKorpi et al.l l|l999l ) have 
included in their 3D numerical simulations the effects of the 
differential rotation and the magnetic field of the galactic 
disk, but again considered a small computational domain 
(500 pc 2 x 1 kpc). In all the simulations above any depen- 
dence of the variables with the Galactocentric distance is 
neglected. The lack of an extended grid makes the entire 
development of a disk-halo-disk cycle very difficult to follow 
because part of the mass is lost through the grid boundaries. 

The simulations bv lFraternali fc Binnevl (|2006h followed 
a different strategy. They considered the whole Galaxy and 
calculated the orbits of gas clouds assumed to be like bullets 
ejected (by SNe II) from the disk that move ballistically up 
into the halo and then return to the disk. However, hydro- 
dynamical effects have been essentially neglected in these 
models. 

In order to study the large scale motion of the ISM and 
the halo gas and follow the complete dynamical and ther- 
mal evolution of the gas lifted by the galactic fountains, the 
whole, differentially rotating Galaxy must be considered. In 
Paper I we have presented 3D hydrodynamical simulations 
of the evolution of a single GF (SGF), i.e. a fountain gen- 
erated by only one OB association with 100 SNe. There, we 
investigated the influence on the SGF dynamics of several 
factors, such as the Galactocentric distance, the Galactic 
rotation, and the presence of a gaseous halo and/or a disk 
neutral hydrogen layer. We have found that a SGF may 
eject material up to ~ 2 kpc before it collapses back in 
the disk, mostly in form of dense, cold clouds and filaments, 



with velocities which are compatible with those of the IVCs. 
Contrary to the common expectation, the gas lifted up by 
the SGF tends to move toward the centre rather than to- 
ward the outskirts of the disk because part of its angular 
momentum is transferred to the halo and/or the upper lay- 
ers of the disk. However, most of the GF gas falls back on 
the disk within a radial distance AR ~ 0.5 kpc from the 
place where the fountain originated. This localized circula- 
tion of disk gas agrees with recent chemical models of the 
Milky Way whi ch assume that the SNII metals enrich only 
the local ISM. (|Cescutti et al.ll2007l ). 

The present paper is devoted to multiple Galactic foun- 
tains (MGFs) - i.e. fountains generated by multiple OB as- 
sociations - in which the gas flow is due to the continu- 
ous star formation in the Galactic disk. Thus, this paper 
represents a natural extension of the analysis of Paper I 
on single fountains. To keep the computational costs un- 
der reasonable limits, we have followed in detail only a lim- 
ited volume of the Galaxy where the MGF occurs; the re- 
maining Galactic volume was mapped at a lower resolution. 
Although a v ery accurate des c ription of the different I SM 
phases (as in Ide Avillel |2000| ; Ide Avillez fc Berrvl l200ll ) is 
hampered by the limited resolution resulting from our ap- 
proach, the global evolution of the fountains can be satisfac- 
torily calculated. We neglected the effect of magnet i c field ; 
as shown by the afo r ementioned papers by Kahn (1 199 lT) ; 
IKahn fc Brettl (|l993l ); Ide Avillez fc BreitschweTdl (|2005l ) it 
does not strongly affect the global fountain dynamics. 

An important motivation to run simulations of MGFs 
is to investigate the mass and energy exchange between the 
star forming disk and the hot gaseous halo surrounding the 
Galaxy. The efficiency of the SN feedback is a key, yet poorly 
understood issue in th eories of galactic star format ion, as we 
will outline later (c.f. iNavarro fc Steinmetd 119971 ). We find 
that a non-negligible fraction of the SN energy is indeed 
injected in the halo (Section 3.4). 

A further reason to consider multiple fountains is to ex- 
plore the interaction with external gas infalling onto the 
disk. Indeed, such gas is needed in current galaxy evo- 
lution models in order to maintain star formation rates 
similar to those observed i n normal spiral galaxies (e.g. 
iTakeuchi fc Hirashital |2000| ; ISemelin fc Combesl |2002| ), as 
well as to explain the obs erved patterns in the chemical evo - 
lution of the Milky Wa y i|Chiappinill2002l ; iGeiss et al.ll2002l ). 
iFraternali et al.l (|2007l ) detected inflowing HI clouds in the 
nearby spiral galaxies NGC 891 and NGC 2403. They ar- 
gue that these clouds are a clue to a substantial and other- 
wise mostly hidden halo gas accretion interacting with GFs. 
While the required accretion rate is comparable to the star 
formatio n rate, the infalling gas mu st have low angular mo- 
mentum (|Fraternali fc Binnevl 120061 ). 

In the following sections, we outline the basic charac- 
teristics and the setup of our numerical model (Sec. 2), then 
examine a reference model for the investigation of multiple 
fountains (MGFs) formation in the Galaxy at a distance of 
8.4 kpc from the galactic centre (Sec. 3) and compare these 
results with a model where the MGFs occur at R — 4.5 kpc 
(Sec. 4). Finally, in Sec. 5, we consider the interaction of 
MGFs with infalling gas from the intergalactic medium and 
in Sec. 6, we draw our conclusions. 
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2 THE MODEL 

The model of the Milky Way, as well as the assumptions rel- 
ative to the SN II explosions belonging to single associations, 
are fully described in Paper I. Here we briefly recapitulate 
them and discuss in detail how we take into account the 
stellar explosions occurring in different OB associations. 



2.1 The Galaxy model 

The ISM in our model is made up of three components - 
namely molecular (Eh), neutral (H i) and ionized (H II) 
hydrogen - a ll follo wing the spatial distribution given by 
IWolfire et all (|2003l ) (see also Paper I). The ISM is ini- 
tially set in rotational equilibrium in the Galactic gravita- 
tional potential due to the summation of the dark matter 
halo, the bulg e and t he disk con t ributi ons, as described by 
iNavarro et al.l (|l996l). iHernquist] (l 99C0. and by a flattened 
King profile (|Brighenti fc Mathews! 1 19961 ) . respectively. A 
hot isothermal gas halo in hydrostatic equilibrium with the 
Galactic potential well is also added. This halo has a cen- 
tral density po,h = 2.17 x 10 -27 g cm -3 and a temperature 
T h = 7 x 10 6 K. 



2.2 Supernovae explosions 

2.2.1 Single associations 

It is known that th e great majority of the massive stars 
form in clusters (e.g. IWard-Thompsonll2002l . and references 
therein). In Paper I, we defined a "single fountain" as a 
fountain powered by the SNe II present in a single clustefl 
The OB association driving the fountain was assumed to be 
totally housed inside a single zone of the highest resolution 
level of our adaptive grid (see section[2]3} • The SNe II do not 
explode all at the same time because of the different mass of 
their progenitors. Thus the energy and the mass delivered by 
all the explosions are released over a time interval of At=30 
Myr (the lifetime of an 8 Mq star, which is the least massive 
SN II progenitor). The stellar explosions inject mass and 
energy at the rates M — TZM^ and L w = IZEq, respectively. 
Here M e j = 16 Mq and Eo = 10 51 erg are the mean mass 
and energy released by a single explosion, respectively; 1Z = 
iVsN/At is the SN II rate, and iVsN is the total number of 
SN II explosions occurring in the association. L w represents 
the mechanical luminosity of the wind driven by the SNe II 
and inflating a superbubble. If the luminosity is larger than 
a critical value L CI , the superbubble breaks out of the disk 
forming a chimney in the ISM and gives rise to the GF (cf. 
Paper I and references therein). In section T2.3I we describe 
as -£/w is implemented into the numerical grid. 



2.2.2 Multiple associations 

At the observed rate (see below), several SGFs may occur 
sufficiently close to each other in time and space to mutually 
interact. This interaction obviously modifies the general dy- 
namics of the SN II ejecta of each single fountain. Following 
with an adequate spatial resolution all the SN II explosions 



For the sake of simplicity, by SNe II we mean here all the core 
collapse explosions. 



occurring in the whole disk is beyond our current compu- 
tational capabilities. Thus, throughout the simulations of 
multiple fountains, we focus in particular on a circular area 
(hereafter referred to as the "active area") of the Galac- 
tic disk of 8 kpc 2 which is resolved up to the maximum 
level of grid refinement (this finest grid extends vertically 
up to z — 3.2 kpc). We allow the SNe II to occur only in 
this area over a period T = 200 Myr. It must be empha- 
sized that the star formation rate (and the associated SN 
rate) is not computed self-consistently with the gas evolu- 
tion. Instead, we assume a mean SN rate adequately scaled 
from the rate 1.4 x 10~ 2 yr _1 relative to the whole Galaxy 
l|Cappellaro et al1ll997t ). It turns out that during the time 
T of our simulations iVtot = 4.2 x 10 4 SNe II are expected to 
explode in the active area of the Galactic disk (see, however, 
the end of this section) . 

Let us assume /(Asn^A^n to be the number of stel- 
lar clusters with a number of supernova progenitors between 
JVsn and Nsn + dNs~N- /(Asn) fol lows a power law distribu- 
tion (|Higdon fe Lingenfelterl 120051 ), 



/(JVsn) « A^ 2 , 



AUn < A^SN < N„ 



We normalize /(Nsn) in order to get 

f JV max 

A'sN/(A r SN)rfA r sN = AW 



(1) 



(2) 



This normalization can be done once N m i n and A^ax have 
been chosen (see below in this section). 
We choose bins in A^sn such that 



N i+1 
SN 



/(AsN)dAsN = 1; 



(3) 



we then assume that such a cluster contains 0.5(ATg N +A r gJ 1 ) 
SNe II. Clearly, for i = 1 we have = A m i r, . 

Once obtained the number and the richness of the clus- 
ters, we place them within the active area. This is made as- 
sociating randomly to each i-th cluster a position Pj within 
the active area and a time tj in the range < tj < T. As the 
time proceeds, the supernovae belonging to the i-th cluster 
start to explode at Pj at the time tj (and continue to ex- 
plode for 30 Myr, as outlined in section [2. 2. 1[) . The sequel of 
the SN II "bursts" is shown in Fig. [T] The energy provided 
by the i-th burst is EoNg N . 

To choose Amin and Af max , we simp ly assume 
A max =2000 from iHigdon fc Lingenfelterl (|2005l ). With re- 
gard to A^min, although in principle it can be as low as 
Amm = 1, we assume N m - m = 30. At a distance of R — 8.5 
kpc, which is the Galactocentric distance of the centre of 
the active area in our reference model (cf. section [3J, GFs 
can form only when L w > L CI = 3 x 10 37 erg s _1 (cf. Paper 
I), that is for clusters richer than A^sn > 30. In order to 
capture the physics of interest here, we have thus adopted 
Nmi n = 30. As a caveat, we stress that with this choice of 
A'min the actual number of SN II explosions occurring after 
T — 200 Myr is nearly Atot/2, i.e., half of those expected. 
Although this is not important from a dynamical point of 
view, it influences the amount of metals set in circulation by 
the fountains. We have to remember this when discussing the 
chemical pollution of the disk. 
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Figure 1. Temporal histogram of the sequence of stellar clusters 
fuelling our models of multiple Galactic fountains. The y-axis 
reports the number of SNe II of each cluster; the x-axis shows the 
time in Myr. 



2.3 Numerical Setup 

We use a modified version of the numerical adaptive 
mesh refinement (AMR ) YGUAZUa hydrodynamical code 
l|Raga et alj|200d . l2002h . The radiative losses are calculated 
assuming solar abundances and by using non-equilibrium 
ionization species abundances (for HI, HII, CII, CIII, CIV, 
OI, Oil and OIII) when T > 10 6 K. For higher tempera- 
tures, where the collisional ioniziation equilibrium is often a 
good appr oximation, we use an equilibr ium cooling function 
taken from lSutherland fe Dop ita (1993). More details can be 
found in Paper I. We enforce the maximum grid resolution 
only within the active volume defined in the previous section 
(given the symmetry of the problem, our computational vol- 
ume encompasses only the "upper" half of the Galaxy). As 
a reasonable compromise between the need of a large active 
surface and a good spatial resolution within it, we chose an 
area of 8 kpc 2 and a minimum mesh size of 50 pc. The other 
three coarser grid levels have mesh sizes of 100 pc, 200 pc 
and 400 pc, respectively. 

The rate of mass and energy injection by SNe II be- 
longing to every single association is described in section 
12.2.11 Because the numerical grid is at rest in space while 
the Galaxy rotates, each OB association occurring in the 
active area moves along a circle on the z = plane. At each 
time we place the energy and mass sources in the numeri- 
cal cell of the finest grid transited by the association at that 
time. We stress that while we assume that all the OB associ- 
ation are located at z — 0, SNII progenitors are observed to 
have a vertical distribution with scale height of ~ 100 pc. As 
our finest grid has a size of 50 pc, in order to investigate the 
effect of the vertical position on the fountain evolution we 
must locate OB associations only two zones above the equa- 



torial plane. Such a small displacement is not expected to 
generate appreciable differences in the calculated flow evo- 
lution. 



3 THE REFERENCE MODEL 

In this section we discuss the reference model (RM) , in which 
the active area is centred at the Galactocentric distance R = 
8.4 kpc. 



3.1 ISM in the disk 

Figure [2] shows the time evolution of the gas column den- 
sity of the ISM of the disk along the direction perpendicular 
to the Galactic plane (the z-direction) . The two dashed cir- 
cles in the panels delimit the region within which the active 
area evolves as the Galaxy rotates. There are several cavities 
whose lifetime is given by the duration of the SN activity in 
each association (~ 30 Myr) added to the time needed for 
the replenishment to occur at the end of such activity. The 



replenishment time is approximately given by t T , 



where r c is the maximum radius attained by the cavities, and 
c is the sound speed of the ambient medium. It turns out 
that r c ~ 0.7 kpc for the largest cavities, while c~2x 10 6 



cm s , and thus t r 



37 Myr. In conclusion, the lifetime 



of the holes into the disk ISM is < 70 Myr. This lifetime and 
the occurrence rate of the OB associations both conspire to 
avoid a substantial overlapping of the holes and the area of 
interest is never totally deprived of its ISM. Given the dif- 
ferential rotation of the Galaxy, the holes become elongated 
with positive pitch angles (angle between the major axis of 
the hole and the tangent to a circle at that galact ic radius), 
as described by the models of lPalous et al, (1990). As a fur- 
ther effect of the disk differential rotation, the holes created 
at earlier times lag behind the recent ones, tending to form 
a sort of queue. This latter effect, however, arises because 
our active region is limited. Had the active area covered the 
whole disk, its appear ance would be similar to NGC 2403 
(|Fraternali et al.ll2002|) and other typ i cal sp irals. 

In a recent paper Boomsma et al.l (|200Sl ) presented a de- 
tailed study of the H I holes and HVCs in the spiral galaxy 
NGC 69 456. It is interest i ng to compare our results to their 
findings. iBoomsma et alj (|2008l ) evaluate an average radius 
of the holes of 0.6 kpc and a maximum lifetime of about 
80 Myr; both these figures are very close to the values ob- 
tained in our simulation. Those authors also find that the 
hole covering factor drops sharply toward the smallest galac- 
tic radii, an effect ascribed to the drop of the H I column 
density as well as to the stronger shear which shortens the 
hole lifetimes. Smaller dimensions and shorter lifetimes at 
smaller radii are also met in our models, and we find that 
the reason is given by the larger pressure of the ambient gas 
which contrasts the bubble expansi on and favours t h e hole 
replenishment (cf. sect . [4j . Finally. IBoomsma et al.l (|2008l ) 
estimate an upper limit of about 10 7 Mq for the average H I 
mass missing from each hole. In our simulation we have 3-5 
holes present at every time (cf. Fig. [2}, while the disk gas 
lifted at any time is a few 10 7 Mq (cf. sect. I3.4 |l. thus giving 
roughly 5 x 10 6 M q per hole, consi s tently with the upper 
limit estimated by iBoomsma et al.l (|2008l ). In conclusion, 
despite our one-fluid approximation of the multi-phase ISM 
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t = 100 Myr 




Figure 2. Face-on view of the active area for the reference model (RM). Column density along the z direction of the ISM of the Galactic 
disk at several times. The two superimposed dashed circles have radii of 6.8 kpc and 10 kpc and they delimit the ring within which the 
active disk area moves. The logarithmic column density scale is given in cm -2 and the spatial scale along the x and y directions is in 
kpc. 



and the galactic model tailored on the Milky Way rather 
than on NGC 69456, our results are in encouraging agree- 
ment with the observational data. 



3.2 Gas circulation 

Figure|3]represents a sort of tomography at different times of 
the column density of the SN II ejecta along the z direction. 
The total column densities are shown in the right panels; 
the left and right panels illustrate the contributions of ejecta 
located below and above z = 800 pc (the height at which the 
disk/halo transition is located for a radial distance of R = 
8.4 kpc), respectively. The dense spots present in the first 
column of the panels are due to recent SN explosions (see, for 



example, the top-left panel). Instead, the more structured 
distribution visible in the middle column is due to ejecta 
trapped and condensed within clouds forming at larger z 
from the lifted disk gas which suffers thermal instability and 
cools down to T = 10 4 K (the minimum temperature allowed 
in our simulations) . These clouds are characterized by large 
column densities, and are visible in Fig. [3] as red spots. 

The cold, dense gas must eventually fall back onto the 
disk. Figure U shows the column density of the clouds with 
negative vertical velocity (the cloud gas is selected as gas 
with T ^ 5 x 10 4 K). This figure clearly demonstrates how 
the gas lifted by the fountains falls back remaining mostly 
within the ring of the disk in which the active area is lo- 
cated. This result confirms what was found in Paper I and 
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Figure 3. Column density along the z direction of the SN II ejecta for the reference model. Each row corresponds to a different time. 
The left and central panels correspond to the column densities calculated in the intervals z < 0.8 kpc and z > 0.8 kpc, respectively, 
where z = 0.8 kpc corresponds to the location of the disk/halo transition. The right column shows the total column densities. The two 
circles have the same meaning as in Fig. [2] The logarithmic column density scale is given in cm -2 and the spatial scale along the x and 
y directions is in kpc. 



indicates that there is no substantial radial dispersal of the 
metals from the Galactocentric radius where they are pro- 
duced. A sim i lar re sult has been also recently reported by 
ISpitoni et all <|2008h . 

Moreover, Figures [3] and [4] show that the small fraction 
of the rising gas escaping to the coronal sector moves radially 
preferentially outward rather than inward. This is due to the 
decrease of the centripetal component of the gravity with 
height and to the conservation of th e angular momentum 
i|Shapiro fc FielJi976tferegmanlll98Ch . The tendency of the 
MGF gas to move outward is however, contrasted by its 
interaction with the hot halo (which is assumed at rest) . The 



fountain transfers part of its angular momentum to the halo 
and thus tends to reduce its radial displacement outward 
(see below). 

A more quantitative insight of the evolution of the gas 
circulation is given by Fig. [5] The top panel shows the loss 
of angular momentum of the SN II ejecta due to the inter- 
action with the halo gas 0. At any time, the angular mo- 
mentum displayed in the figure is normalized to the total 



2 We point out that the time evolution of the relative angular 
momentum of all the gas lifted by the fountains is quite similar 
to that relative to the SN II ejecta alone. 




Figure 4. Column density along the z direction of the gas falling back to the disk, for the reference model. The times are the same as in 
Fig. [21 as well as the meaning of the two circles. As apparent from the figure, we consider only gas condensed into clouds and filaments. 
It is also clear that the clouds remain essentially within the circular sector in which the active area is located. The logarithmic column 
density scale is given in cm' 2 and the spatial scale along the x and y directions is in kpc. 



angular momentum injected by the SN II up to that time. 
The middle panel illustrates the fraction of the SN II ejecta 
located at different heights, while the bottom panel shows 
the fraction of ejecta moving outside of the boundary of the 
active area and testifies that the amount of gas moving out- 
ward is larger than that moving inward. The latter starts 
to increase after ~ 70 Myr, when a substantial amount of 
the lifted gas starts to fall back (see section [3T4J ; having lost 
part of its angular momentum, this gas tends to return at 
a Galacocentric radius smaller than that at which it started 
its journey. 

It is interesting to compare Fig. [S] with the analogous 
figure in Paper I (Fig. 5) in order to highlight similarities and 
differences with the SGF case. The relative loss of the an- 



gular momentum of the fountain gas is essentially the same. 
This must be expected as in MGFs a larger amount of foun- 
tain gas interacts with a larger amount of halo gas, leaving 
the loss of specific angular momentum essentially unaffected 
compared to SGFs. For analogous reasons, also the rela- 
tive amount of gas lifted above different heights is similar; a 
larger amount of energy, in fact, moves a larger amount of 
gas to similar altitudes. Obviously, as the SN explosions keep 
going for a longer time, new gas persists at larger heights, 
at variance with the SGF case. Finally, we stress that in 
the case of SGFs on the long term the fraction of gas mov- 
ing outward is smaller than that moving inward, at variance 
with the MGF case. However, a glance to Fig. 5 of Paper I 
reveals that, initially, also in the SGF case more gas moves 
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outward rather than inward. This different behaviour is due 
to the fact that in the present model the MGF flow is con- 
tinuously supplied, while in the SGFs the SN engine stops 
its activity after 30 Myr and the fountain gas, being no more 
replenished with fresh angular momentum, is facilitated in 
its centripetal displacement. 

In any case, it is quite clear that most of the gas remains 
essentially at the same Galactocentric distance where it has 
been created. 



3.3 Cloud properties 

We now discuss in more detail the structure of the clouds 
forming from the gas powered by the MGFs. In particular, 
the ISM lifted by the fountains condenses essentially at the 
top of its trajectory, and than falls back with velocities < 100 
km s _1 . ~ 130 Myr after the beginning of the simulation an 
equilibrium is established between ascending and descend- 
ing gas (cf. sect. !3.4) l. After this time, the gas lifted above 1 
kpc has a mass 1.5 x 10 7 Mg, 75% of which is condensed in 
dense filaments cooled down to T = 10 4 K (see Fig.0. The 
clouds are composed essentially of gas originally from the 
disk (~96%), and only a slight contribution (4%) from the 
halo gas that was compressed and cooled after interacting 
with the fountains. As we have found in Paper I, the SN 
ejecta plays a negligible role in the mass budget of the foun- 
tains. After 130 Myr, 1.3 x 10 4 SNe II have exploded and 
delivered 2.1 x 10 Mq of ejecta, but only half of it enters 
within the computational domain (the rest is injected in the 
"lower" half of the halo, which was not considered in our 
simulations, for symmetry reasons, as remarked in Sec. 2). 
Nearly 50% of the ejecta remains trapped within the disk, 
35% is encompassed within the evolving clouds, and 15% 
floats over the disk as hot, diffuse gas. Thus the SN ejecta 
contribute for less than 0.4% to the clouds mass. For this 
reason, the chemical composition of the clouds is unaffected 
by the metals produced by the SNe II driving the fountains. 
As each supernova d elivers on average 3 Mq of metals (cf. 
iMarcolini et al.ll2007h . a total mass of 2 x 10 4 M of heavy 
elements is ejected by the fountain within the computational 
domain. As a result, the metallicity increment in the clouds 
due to the freshly delivered metals is AZ = 5.7 x 10 -4 , 
rather negligible compared to the solar metallicity of the 
ISM. This conclusion holds also if one considers the contri- 
bution of SNe II belonging to OB associations hosting less 
than 30 SN progenitors, and representing half of the total 
number of supernovae (see section I2.2.2|l . 

In conclusion, the clouds forming in our simulation have 
solar metallicities and velocities lower than 100 km s -1 , and 
are th erefore likely to be associated to IVCs rather than 
HVCs (|Wakker et "alj|2008h . This result is in agreement with 
that found in Paper I when considering single fountains. 



Figure 5. Illustration of several quantities for the reference 
model. Top panel: angular momentum evolution of the SNII ejecta 
normalized to the angular momentum injected by the SN II ejecta. 
Middle panel: fraction of the mass of SN II ejecta located at 



— 3.4 Disk-halo energy exchange 



Other details of the gas flows of the reference model are given 
in Fig. [7] and are intended to shed light on the feedback 
between the disk and the hot halo. This energy exchange 
regulates the evolution of the halo gas and influences the 
star formation history in the disk. It is also relevant in the 
related issue of the so-called "over-cooling" problem in cos- 
mological simulations of galaxy formation, which generally 
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overestimate the amount of radiatively cooled gas in galax- 
ies. Energy feedback from SNe and/or AGN, albeit poorly 
understood, is thought to be a crucial ingredient to reconcile 
hierarchic al models of galaxy formation with the observa - 
tions (e.g. iNavarro fc Steinmet3ll997t ICroton et al]|200rj) . 
Supernov ae alone are probably incapable to fully solve the 
problem (jTornatore et al.l 120041 ). but is nevertheless impor- 
tant to quantify the SN energy which heats the hot halo gas 
to correctly estimate its cooling rate and the disk-halo mass 
exchange. 

As in our model the amount of thermal energy of all 
the halo is much larger than the energy injected by the SNe, 
a direct evaluation of its variation due to the stellar explo- 
sions would be rather inaccurate (also because of the per- 
turbations generated at the grid boundaries, as discussed in 
Paper I). We thereby check the energy balance of the halo- 
fountains interaction as follows. In the top panel of Fig. [7J 
we report the thermal energy of the disk gas that is lifted by 
the fountains and not condensed into clouds. This rarefied 
gas endures over the disk and its thermal energy is incorpo- 
rated in the hot halo. On the other hand, some coronal gas 
cools and condenses over the clouds created by the MGFs. 
The thermal energy radiated by this cold gas is also shown 
in the top panel of Fig. [7] At the end of the simulation, the 
halo has gained an energy twice larger than that radiated 
due to the interaction with the fountains, with a net increase 
of thermal energy of ~ 10% of the energy released by all the 
SNe. This value for the SNe heating efficiency is expected 
to remain stable once an approximately steady regime is 
established between rising and falling gas. 

The achievement of this regime is visible in the middle 
panel of Fig. [JJ that shows the fraction of rising (solid line) 
and descending (dashed line) fountain gas. At early times 
only the rising gas is present. After nearly 40 Myr the first 
clouds form and fall back. A dynamical equilibrium between 
ascending and descending gas is established after 130 Myr. 
The same result is illustrated in the bottom panel in terms 
of the absolute values of the rising and falling masses. 

As pointed out, we cannot make an exact evaluation 
of the energy balance between the fountains and the coro- 
nal gas. In the analysis above (see the top panel of Fig. 7), 
we neglected the gain of the kinetic energy of the hot gas 
interacting with the galactic fountains and the change in 
the potential energy of the whole halo. Also, we neglected 
the radiative losses of the hot gas, possibly enhanced by the 
compression due to the fountain expansion. However, the 
flow of the halo gas remains subsonic and its cooling time 
is longer than the time of the simulation. We thus conclude 
that a net energy gain is attained by the halo during its 
interaction with the MGFs. 



4 MODEL AT R=4.5 KPC 

In this section we explore the influence of the Galactocentric 
distance on the gasdynamics of the MGFs. In particular, we 
have run a model similar to the RM, but with the active 
area centred at R = 4.5 kpc. Here the gas of the disk is 
denser and has a larger effective vertical scale heigh H a s . As 
a consequence, the disk-halo transition occurs at z = 1.35 
kpc, instead of z = 0.8 kpc as in the RM (see also Paper I). 
The larger density of the ambient medium makes more 




o 



Figure 6. Distribution of the disk gas at t = 130 Myr at the 
plane orthogonal to the Galactic disk and passing through the 
active area for the reference model. Isodensity curves of the ejecta 
are over-imposed, highlighting the fountain pattern and the cloud 
formation. Since the density of the ejecta is much smaller than 
that of the ISM, it is represented in a different scale in order to 
be visible in the figure. The logarithmic density scale of the disk 
gas is given in g cm -3 and the spatial scale along the x and z 
directions is in kpc. 



difficult the expansion of the SN bubbles through the plane 
and accelerates their backfilling. The reduced size and life- 
time of the holes also reduce their probability of merging. 
This is apparent in Fig. [8] where the number and the size 
of the holes are smaller than in the reference model (cf. Fig. 

rj. 

In Paper I it was shown how the vertical distribution 
of the disk gas influences the dynamics of a SGF and the 
shape of the chimney carved by the fountain. In particular, 
a single fountain located at R = 4.5 kpc gives rise to a sort of 
well collimated outflow that crosses the disk/halo transition 
only barely. On the contrary, given the lower density and 
scale height of the local gas, the same fountain located at 
R = 8.5 kpc experiences a well defined break out. In the case 
of MGFs these differences are greatly reduced because the 
energy powering the gas flow in this case is much larger than 
the critical luminosity L CI at any Galactocentric distance. 

A somewhat larger alignment is still present at R — 
4.5 kpc also in the case of MGFs, as apparent when one 
compares Fig. [3] with Fig. [9] (where the different position of 
the disk/halo transition has been taken into account). The 
amount of diffuse gas spread outside the circular sector is 
smaller than in the RM. The same conclusion can be reached 
looking at Fig. 1101 The falling gas is clearly less radially 
dispersed compared to the reference model (cf. Fig. 3| . Note 
that, since in the present model the circumference covered 
by the active area is shorter, after 200 Myr a complete ring 
is described in all the panels of the bottom raw of Fig. [TUJ 

As for the RM, more quantitative information can be 
obtained looking at the Fig. 1111 The loss of angular mo- 
mentum is essentially the same as in the RM (top panel). 
Comparing to the SGF evolution (Fig. 9 in Paper I), we re- 
mark that in this latter case the angular momentum, after 
an initial drop, slowly goes up again. As discussed in Paper 
I, this is due to the fact that a larger fraction of the SGF gas 
falls back onto the thicker disk gaining new angular momen- 
tum. This effect is not present in the MGF model because, 
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t = 50 Myr 



t = 100 Myr 




Figure 8. The same as in Fig. [2] but for the Galactocentric distance R = 4.5 kpc. 



as stressed above, given the larger energy of the fountain, 
the disk thickness does not play a significant role 

The vertical extent covered by the gas (middle panel 
of Fig. HHI is somewhat lower than in the RM as the den- 
sity and the scale height of the disk are larger in this case. 
The differences with the SGF case are analogous to those 
discussed for the RM, and do not repeat them here. 

The smaller radial dispersion of the MGF gas empha- 
sized above is illustrated in the bottom panel of Fig. [11] 
showing that the amount of gas moving radially outside the 
active area in nearly halved compared to the RM. 

Finally, from the top panel of Fig. [12] we see that in 
this model the halo gas acquires 10% of the SN energy at 
the end of the simulation, i.e., the same amount as in the 
RM. Also the time at which the quantity of ascending and 



descending gas become comparable is similar to that of the 
RM (middle and bottom panels of Fig. I12[l . 



5 INTERACTION WITH ACCRETING GAS 

As pointed out in Sec. 1, the gas flow of the fountains is likely 
to sweep up cold gas, but it is unclear whether this gas arises 
from cooling of the hot halo or from co l d stre ams falling 
from large distances. iFraternali fc Binnevl (|2008h argue that 
this latter possibility is the most realistic. If the cold falling 
gas is clumpy, its clouds cannot be more massive than 10 7 
Mq , otherwise th e solar neighbourhood would be heated 
by th eir passage l|Lacev fc Ostrikerl 1 19851 ; iToth fc Ostrikerl 
1992). Actually, while HVCs sur rounding our Galax y have 
masses of order 10 6 " 7 M (e.g. IWakker et alj|2008t ). deep 
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Figure 9. The same as in Fig. [3] but for the Galactocentric distance R = 4.5 kpc. In this case the left and central panels give the 
column densities calculated in the intervals z < 1.3 kpc and z > 1.3 kpc, respectively, since at this Galactocentric distance the disk/halo 
transition occurs at z = 1.3 kpc. 



observations around M31 reveal th e presence of numerou s 
HVCs with masses down to 10 4 M Q (jWestmeier et alj|2005l) . 

The study of the possible accretion of gas onto galax- 
ies is important for several reasons. It is generally be- 
lieved that accretion must replace the gas consumed by 
star formation which has been remarkably constant in 
the solar neighbo urhood over the Galaxy lifetime (e.g. 
iBinnev et afll200Ch . Another reason for studying the galac- 
tic accretion is linked to the mounting evidence of massive 
haloes of neu tral and ionized gas surrounding nearby spiral 
galaxies (e.g. ISancisi et aT1l2008l ). These thick disks rotate 
more slowly than th e thin disks and show inflow motions. 
iFraternali fc Binnevl ()2008h have shown that these haloes 
cannot be sustained only by galactic fountains. According 



to them, the gas from the fountains interacting with a pre- 
existing hot corona would make it to co-rotate with the disk 
after a very short time (< 1 Myr) which is not consistent 
with the observations. This result led them to conclude that 
a substantial accretion of low angular momentum material 
from the intergalactic medium would be required in order 
to assure a slower rotating halo. 



In the model of IFraternali fc Binnevl (|200&t ) the clouds 
are treated as bullets moving ballistically, and hydrodynam- 
ical effects are absent or only roughly represented. In order 
to investigate the issues above more realistically from an 
hydrodynamical point of view, we have simulated the in- 
teraction of the MGFs with infalling gas considering two 
different models, one in which the gas infall is described as 
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Figure 10. The same as in Fig. [4] but for the galactocentric distance R = 4.5 kpc 



a continuous drizzle, and another in which a single, denser 
cloud is assumed to fall toward the active area. 



5.1 The "drizzle" model 

This model is similar to the reference model, but we assume 
that a gas of density n = 10~ 3 cm" 3 drops from the "top" 
boundary of the computational grid with a vertical veloc- 
ity u z = —150 km s _1 . Such a flux extended over all the 
Galactic disk would correspond to an accretion rate of 5.2 
Mq yr -1 (from both sides). This is much higher than the 
accretion rate M ~ 0.2 Mq yr" 1 associated with the HVCs, 
and also higher than the rate of ~ 2 Mq yr~ needed to 
sustain the star formation in the solar neighbourhood. How- 
ever, the accretion rate we assume is comparable to that in 
iFraternali fc Binnevl i|200Sf ) (which was adopted for the case 



of NGC 891) and this will allow a direct comparison with 
their work. 

As we consider the fountain activity only in a limited 
region of the disk plane, we let the rain fall from the top 
boundary only over the ring sector on the disk where the 
active area moves (cf. Fig. [2} • We have considered only this 
region because most of the gas set in motion by the fountains 
remains within this sector (as we have seen in section [3}, so 
that all the features of the interaction of the fountains with 
the rain can be captured. 

As the infalling gas interacts with the hot halo and with 
the ascending gas of the fountains, vortices of ~ 2.5 kpc 
form and the density of the gas above the disk fluctuates, as 
shown in Fig. 1131 These vortices are likely an artifact of the 
method we use to inject the accreting gas. However, these 
features do not alter the flow generated by the fountains, 
being the vortices energetically negligible compared to the 
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Figure 7. Illustration of several quantities of the reference model. 
Top panel: amount of thermal energy of the disk gas transferred 
to the gaseous halo (solid line), thermal energy lost by the halo 
gas condensed onto the clouds (dashed line); see text for details. 
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Figure 11. The same as in Fig. [5] but for the model with the 
active area centred at the Galactocentric radius R = 4.5 kpc. 
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Figure 12. the same as in Fig. [7] but for the model with the 
active area centred at the Galactocentric radius R = 4 kpc. 




Figure 13. The same as in Fig. [6] but for a model considering 
a continuous external gas infall over the active area (the drizzle 
model). 



energy injected by the SNe. Figure [13] corresponds to the 
same time and grid slice of the snapshot depicted in the RM 
of Fig. EI Comparing both models, we see that in the drizzle 
model the disk gas set in motion by the fountains reaches 
lower heights and tends to lag behind the active area (toward 
the right hand side in the figure). 

A more quantitative analysis is given in Fig. [14] The 
angular momentum loss (top panel) is more than doubled 
relative to the RM. This explains the larger lag of the foun- 
tain flow relative to the disk, as noted above. The middle 
panel quantifies the vertical distribution of the fountain gas; 
it is essentially similar to the reference model up to z = 1.2 
kpc (cf. Fig. [5]l, but less gas is found above z — 2.4 kpc. 

Once an overall balance is reached between ascending 
and descending gas, as in the RM, 25% of the lifted gas 
remains rather diffuse, while 75% of it condenses in dense 
filaments. The total mass of these clouds is 7 x 10 6 M©, 
i.e., nearly 60% of what was found in the RM; such reduc- 
tion is related to the lower vertical excursion of the clouds 
and their shorter "lifetime". It is interesting to analyze the 
"composition" of the clouds. As in the RM, the contribu- 
tion to their mass by the SN ejecta amounts to a negligible 
0.3%. This time, however, the amount of disk gas is lower 
(87.9%), while the fraction of "external" gas is three times 
larger (11.8%). This result is obviously due to the addition 
of the rain material to the halo gas. This denser ambient 
medium cools more easily when it interacts with the clouds 
and condenses onto them more efficiently. 

The accretion ma ss onto the clouds is a cru cial ingre- 
dient in the model of iFraternali fc Binnevl (|2008l ). In their 
model, they integrated the orbits of the clouds that are 
ejected from the disk of a spiral galaxy and move up ballis- 
tically to the halo and then return to the disk. As the clouds 
accrete mass from external infalling gas, their velocity de- 
creases due to momentum conservation. These authors have 
then shown that, if the thick disk of neutral gas of the spiral 
galaxies is made up of clouds lifted by galactic fountains, 
the accretion of external gas with low angular momentum 
could explain why the observed rotational curve of these 
disk galaxies decreases with z. It is interesting to compar e 
our results with those found by IFraternali fc Binnevl (2008). 

In order to evaluate the mass accretion onto our clouds, 
let us define mo as the "initial mass of the clouds" (i.e., the 
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mass of the disk gas lifted by the fountains that has con- 
densed into clouds), and 5m as the amount of ambient gas 
(i.e. the gas of the halo and the drizzle) that is swept by the 
clouds and increase their mass up to m — moSm. From our 
simulations, we obtain Sm/mo = 0.13; this result is in rea- 
sonab le agreement with what found by IFraternal i fc Binnevl 
|200Sl ) in the case of NGC 2403 (cf. their Fig. 6). 

We can pursuit further the analysis above comput- 
ing the specific accretion rate a = rh/m. To this end we 
must evaluate the clouds' lifetime tn in order to work out 
q = 5m/(mot\{). From Fig. [TS]we note that ascending and 
descending gas reach a dynamical equilibrium after t > 70 
Myr. This is thus the characteristic time spent by a cloud 
during its journey, we obtain a ~ 1.9 Gyr , which is in 
good agreement with the value a ~ 2 Gyr -1 found by 
iFraternali fc Binnevl l|2008l ). Although this agreement must 
not be super-estimated given the approximations involved 
in the previous calculations, yet we consider it significant. 

We now evaluate how the external mass accretion onto 
the clouds affects their rotational velocity. Obviously, given 
the limited extension of the active area in our model, we 
cannot build a thick disk of neutral hydrogen out of it, nor 
reproduce the entire disk rotational curve. However, some 
interesting conclusions can be drawn from the simulations 
above. The time evolution of the rotational velocities at z = 
1.8 kpc are illustrated in the top panel of Fig. [TBJfor both 
the reference and the drizzle models. They give the mass 
weighed mean velocity of a "slice" of clouds at z = 1.8 kpc. 
Likewise the angular momentum (cf. the top panels of Figs. 
[5l and [T4)) . the values of the velocities tend to stabilize after 
an initial drop. 

It is an easy matter to verify whether the difference in 
the cloud rotation velocity between the drizzle model and 
the RM is entirely due to the difference in the amount of 
mass increment (as suggested by the ballistic model of Fra- 
ternalli & Binney 2008), or not. In that case we would have 

/ 1 + §2H 

^rot _ mo / .\ 

Vrot 1 + ' 

where the prime quantities refer to the drizzle model. From 
Fig.[l6lwe find that at t = 200 Myr, v' TOt /v Iot = 0.83. On the 
other hand, computing the amount of accreted matter onto 
the clouds from the simulations we find that the right hand 
side of equation (4) gives a larger ratio v' Iot /v TO t = 0.92. 
This means that only part of the rotation velocity drop is 
due to the accretion of external material to the clouds. A 
rough estimate of the total momentum flux of the system 
(that includes both the clouds momentum flux and the halo 
gas and ram pressures) indicates that the remaining drop of 
the rotation velocity of the clouds is due to their interaction 
with the hot halo which in the case of the drizzle model is 
much denser due to the external accretion. 

In conclusion, though our simulations are restricted to a 
small area of the Galactic disk, and therefore are not suitable 
to reproduce the rotational curve of the entire thick disk or 
of the hot halo, they seem to indicate that the presence of 
an external gas infall may help to slow down the rotation of 
the gas in the clouds (see Fig. 16) and thus the amount of 
angular momen tum that is transferred to the coronal gas, 
as suggested by IFraternali fc Binnevl i|2008f ). 

A final remark is in order with regard to the discus- 



sion above. In the model of IFraternali fc Bin ney (2008]) the 
clouds are described as pellets that preserve their identity 
during all their journey until the return to the disk. Let us 
define x as the ratio between the density of the cloud and 
that of the hot halo, and the flow time as ifl ow = D/v, where 
D and v are the cloud dia meter and velocity, respectively. 
IFraternali fc Binnevl (|2008l ) estimate that typically x ~ 300, 
D = 100 pc and v = 100 km s" 1 . With these fi gures, it turns 
out that the orbital time of a cloud is of the same order of 
the drag time idrag = X^flow ~ 300 Myr. However, a cloud 
can be destroyed by Kelvin-Helmholtz (K-H) instabilities in 
a time scale £kh = Y° 5 t fiow ~ 17 Myr. Thus, in the model 
of IFraternali fc Binnevl (|2008l ) the cloud material is likely 
to linger above the disk as diluted gas. In our models, in- 
stead, the clouds never reach very large densities, and for the 
densest structures we obtain x ~ 100 and D ~ 500 pc. As a 
consequence, it is £kh ~ 50 Myr, which is of the same order 
of the cloud lifetime. Simulations with higher spatial resolu- 
tion would be required in order to follow appropriately how 
the cloud evolution is influenced by K-H instabilities. Such 
more demanding simulations are planned in a near future. 

5.2 The "cloud" model 

In this section we study the effect of an external cloud falling 
toward the galaxy on the fountain gas dynamics. We do not 
follow the sub sequent interact ion between the cloud and the 
disk gas (see iBregmanl |2004| . and references therein for a 
review on this subject). We consider a cloud entering the grid 
with a velocity « z = —150 km s _1 from the top boundary, 
with mass M = 3.5 x 10 6 Mq, which is consi stent with 
the estimate for Complex C (|Miller et all 120081 ) . We also 
assume T = 10 4 K and density n = 0.2 cm~ 3 . The cloud 
is assumed to be initially cylindrically shaped with both 
diameter and height equal to 800 pc. Figure [T7] shows a 
sequel of snapshots taken at different times and illustrating 
the dynamical interaction of the falling cloud with the rising 
gas of the MGF. 

The figure shows that the cloud actually interacts only 
with a fraction of all the MGF gas. In fact, the cloud diam- 
eter is ~ 1.2 kpc (the cloud expands somewhat during its 
trajectory before entering in equilibrium with the ambient 
gas), and its covering factor relative to the active area is 
~ 0.15. As a consequence, the general characteristics of the 
fountains are not greatly altered by the cloud interaction, as 
clearly illustrated in Fig. 1181 



6 CONCLUSIONS 

In Paper I we studied the evolution of a single fountain, i.e. 
the flow generated by a single SN II association, in order 
to understand the influence of a number of factors (such as 
the Galactic rotation, the distance from the centre of the 
Galaxy, and the presence of a hot Galactic halo) on the 
fountain gas circulation and on the properties of the clouds 
forming in this process. In the present paper, we have instead 
considered multiple GFs, i.e., the action (and interaction) of 
several GFs powered by SN II associations sufficiently close 
one to another in space and time. As in Paper I, we focused 
on a model located at R = 8.4 kpc, the solar Galactocentric 
distance. Our findings can be summarized as follows: 



16 C. Melioli, F. Brighenti, A. D'Ercole, E.M. de Gouveia Dal Pino 



24.4 




Figure 17. The same as in Fig. \E\ but for the cloud model, where an external cloud falls over the MGFs. The four panels represent 
snapshots taken at t = 130 Myr (top-left), t = 136.5 Myr (top-right), t = 143 Myr (bottom- left), t = 149.5 Myr (bottom-right). As the 
Galaxy rotates, the falling cloud appears to move rightward. 



(i) The spatial and temporal random sequel of SNe II 
produces holes on the disk, giving "punctured" look to it, 
qualitativel y similar to what is ob served in many galaxies, as 
NGC 2403 (|Fraternali et alj|2002h . The detailed "topology" 
of the gaseous disk depends on several factors such as the 
SN II rate, the size of the holes (which in turn depends on 
the clustering of the SNe and the local density of the ISM - 
See Section 4), and the speed at which they disappear once 
the SNe responsible of their formation stop exploding. 

(ii) As expected, a larger amount of extraplanar gas can 



be found at larger heights relative to a SGF model. How- 
ever, and somewhat more unexpectedly, the fountain flow 
remains radially localized, just as in the SGF case. The ra- 
dial dispersion of the MGF gas, instead, is smaller than in 
the case of a SGF. Given the absence of significant radial 
flows, the fountains play a negligible role in shaping the ra- 
dial profile of the disk chemical abundance, as required by 
chem ical models which ign ore the hydrodynamics of the ISM 
(e.g., (|Cescutti et alj|200% ). 

(iii) Multiple fountains can form extraplanar gas at larger 
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Figure 14. The same as in Fig. [5] but for the drizzle model. 



Figure 15. The same as in Fig. [7] but for the drizzle model. 
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Figure 16. Time evolution of of the circular velocities of the 
clouds at z = 1.8 kpc for the reference model (thin line) and the 
drizzle model (thick line). 



heights compared to SGFs. However, despite their higher 
energy budget, the MGFs are not able to produce clouds 
at an altitude larger than ~ 3.5 kpc above the Galactic 
plane. Importantly, these clouds are mainly made of disk 
ISM (thus sharing the solar metallicity), rather than cooled 
SN ejecta or halo gas. They fall back toward the Galactic 
plane with velocities lower than 100 km _ . From the con- 
siderations above we conclude that the gas condensations 
rising in our simulations can be identified with IVCs (rather 
tha n HVCs). In this sce nario, the IVCs have a local origin 
(c.f. IWakker et aLll200Sl ). 

(iv) After the break out the multiple fountains protrude 
into the hot gas surrounding the Galaxy and an energy ex- 
change is established between them and the hot gas. The 
diffuse fountain gas which does not participate of the cloud 
formation stays over the disk and becomes part of the halo 
gas. Its thermal energy remains in the hot halo. On the 
other hand, part of this coronal gas cools radiatively and 
condenses on the clouds. We find that the final balance is 
in favor of the halo, which is heated by the SNII with an 
efficiency of ~ 10%. While this moderate value implies that 
SN feedback is probably not able to solve the "overcooling" 
problem that arises in the cosmological context of galaxy 
form ation (e.g. iNavarro fc Steinmet j[l997l ; iTornatore et alj 
120041 ) it is more that enough to mantain an extended hot 
gaseous halo and cold/ warm extraplanar gas around typical 
spiral galaxies. 

(v) Our hydrodynamical simulations indicate that the 
extraplanar gas generated by the fountains accretes mass 
from the infalling IGM at a rate similar to that inferred by 
iFraternali fc Binnevl |200Sl) . Though our simulations are re- 
stricted to a small area of the Galactic disk and therefore, 
are not suitable to reproduce the rotational curve of the en- 
tire thick disk, they show that the presence of an external 
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Figure 18. The same as in Fig. [7] but for the cloud model (thin 
lines) compared with the reference model (thick lines). 
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gas infall may help to slow down the rotation of the gas in 
clouds (Fig. 16) and thus the amount of angular momentum 
that is transferred to the halo. However, higher resolution 
simulations are still needed in order to explore the role of the 
K-H instabilities upon the clouds frag mentation and disinte- 
gratio n which have been neglected by iFraternali fe Binnevl 
(2008) and which are not properly resolved in our simula- 
tions. 

Finally, we considered the accretion of single, isolated 
clouds. Our simulations show that this kind of accretion does 
not alter the general characteristics of the MGFs that arise 
from the disk. 
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